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Entropy of monopoles from percolating cluster in quenched SU(2) QCD* 
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The length distribution and the monopole action of the infrared monopole clusters are studied numerically 
in quenched SU(2) QCD. We determine the effective entropy of the monopole currents which turns out to be 
a descending function of the blocking scale, indicating that the effective degrees of freedom of the extended 
monopoles are getting smaller as the blocking scale increases. 



1. INTRODUCTION 

The dual superconductor picture of the 
QCD vacuum is based on the existence of Abelian 
monopoles which appear naturally in an Abelian 
gauge. There are various indications that the mo- 
nopoles arc responsible for the color confinement 
(for a review, see Ref. |2j). One of the most im- 
portant results is the observation of the monopole 
condensate in the confinement phase |3I4| . 

The monopole trajectories form clusters. Typ- 
ically, each lattice configuration contains many 
finite-sized (ultraviolet) clusters and one large 
percolating (infrared) cluster |5l(il7j . The IR clus- 
ter - which occupies the whole lattice - repre- 
sents the monopole condensate |5]. The tension 
of the confining string gets a dominant contribu- 
tion from the IR cluster 6 . In the deconfmement 
phase the IR cluster disappears |5l(i| . as expected. 

The UV and IR monopole clusters were studied 
previously in Refs. |6I7I8I9I10| . Below we inves- 
tigate the action and length distribution of the 
infrared monopole cluster for various lattice sizes 
and scales b at which the magnetic charge is de- 
fined. Using the action and length distribution 
we investigate the entropy of the IR clusters. 
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2. MODEL 

In our simulations we use the Wilson action, 
S(U) = -f TrlJp. We work in the Maximal 
Abelian (MA) gauge |^. The Abelian gauge 
field, 6, is determined as a phase the diagonal 
component of the SU(2) link variable, p (s) = 
arg l/"(s). The Abelian field strength, 9 pv {s) £ 
(— 47r, 47r), is decomposed into two parts, tlv {s) = 
0^(s) + 27rm M;y (s), where O^s) € [—it, n) is in- 
terpreted as the electromagnetic flux through the 
plaquette and m^ v {s) can be regarded as a num- 
ber of the Dirac strings piercing the plaquette. 

The elementary monopole currents are defined 
as follows |T2], k p (s) = l/2e^ pa d u m pa (s + fi), 
where d is the forward lattice derivative. To study 
the monopole charges at various scales we use the 
n 3 extended monopole construction [3], 

n-l 

k l n) (s) = kn(ns+(n-i)p, + iv + jp + l6-). 

i,j,l=0 

The extended monopole is defined on a sublattice 
with the spacing b = na, where a is the spacing 
of the original lattice. Both elementary and ex- 
tended monopole charges are conserved and quan- 
tized. 

We used the standard Monte-Carlo procedure 
to generate 1000-3000 configurations of the gauge 
field for each value of j3 = 2.1 ~ 2.6. To fix the 
MA gauge we used either the usual iterative algo- 
rithm (for 24 4 and 48 4 lattices) or the Simulated 
Annealing method |13| with five Gribov copies 
(for 32 4 lattice). 
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3. MONOPOLE ACTION 



4. LENGTH DISTRIBUTION 



To get the monopole action we integrate out all 
degrees of freedom but the monopole ones. Fol- 
lowing Ref. [3] we generate the SU(2) configura- 
tions, fix the MA gauge, then get the configura- 
tions of the IR monopole clusters as it was de- 
scribed in the previous section. Then we use the 
inverse Monte-Carlo method to get the monopole 
action. 

The monopole action can be represented |3I15| 
as a sum of the n-point (n > 2) operators Sf. 



fi Si [k] , 



where fa are the coupling constants. We 
adopt only the two-point interactions in the 
monopole action (i.e. interactions of the form 
Si ~ k fJ ,(s)k l j / /(s')). A detailed description of the 
Si[k] interactions can be found in Ref. [16|. 



• 24 , all trajectories 

■ 48 , all clusters 

O 24 4 . IR cluster 

□ 48 4 . IR cluster 
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Figure 1. The self-interaction coupling fi vs. b. 

One can find that the monopole action is pro- 
portional with a good accuracy to the length of 
the monopole loop, 5 mon = fL + const. The dom- 
inant term in the monopole action corresponds 
to the most local self-interaction of the monopole 
currents, Si[k] — J2 S fj.^fii 3 )- We compare the 
parameter /i for the monopole action of the IR 
cluster and for the action associated with the 
whole monopole ensemble in Figure ^ 

One can notice that the coupling constant /i 
shows scaling for large b = na in agreement with 
Ref. 3 . This coupling is independent of the lat- 
tice volume and for large blocking scales b the 
type of the ensemble (the IR cluster or the whole 
ensemble) is not essential for determination of /i . 
However, at small b values, by/cr < 0.5, the type 
of the lattice ensemble becomes important. 



The distribution of the ultraviolet clusters 
was studied both numerically |7I1()| and analyti- 
cally |9I14| . The distribution can be described by 
a power law Djjv oc L~ T , where the power r is 
very close to 3, Ref. jjj. This behaviour indicates 
that the monopoles in UV clusters are randomly 
walking objects |2J- In our simulations we are con- 
centrated on the IR cluster because the IR cluster 
is important for the confinement of quarks. 

The length distribution function, D(L), is pro- 
portional to the weight with which the particular 
trajectory of the length L contributes to the par- 
tition function. From the previous section it is 
clear that monopole action contributes to D(L) in 
a form of an exponential factor, oc e~f L . The en- 
tropy of the monopole trajectory also contributes 
to the monopole length distribution. The entropy 
contribution is proportional to /i L (with n > 0) 
for sufficiently large monopole lengths, L. Thus 
the distribution of the monopole trajectories in 
infinite volume must be described by the function 



DfX(L) oc fi L ■ e~ 



7 = ln/i-/. (1) 



In this equation we neglect a power-law prefactor 
which is essential for the distribution of the ultra- 
violet clusters. In the finite volume there appears 
a finite-volume suppression factor and the total 
distribution can be described as a Gaussian |16| : 

D IR (L) = const. exp{-a(6, V)L 2 + ~/{b)L} . (2) 

The peak of this distribution, L max — "/(b)/ 
2a(b, V), is expected to be proportional to the 
volume of the system, V , to insure the finiteness 
of the IR monopole density, pm — L rnax /V, in 
the thermodynamic limit, V — ► oo. Thus we con- 
clude that a(b, V) = A(b)/V, where A(b) is a cer- 
tain function. Therefore in the thermodynamic 
limit the parameter a vanishes and the distribu- 
tion J5J is reduced to Eq.QJ, as expected. 

To get the monopole entropy we need to know 
the parameter 7, Eq. QJ. We fit the monopole 
loops distributions by the function J5J) and use 
the bootstrap method to estimate the statistical 
errors. We obtain that the parameter 7 — shown 
in Figure |3 — scales with b — n ■ a and is inde- 
pendent on the volume of the system. 
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Figure 2. The fitting parameter 7 as functions b 
for various lattice volumes and blocking steps, n. 

5. MONOPOLE ENTROPY 

The knowledge of the distribution and the 
monopole action allows us to determine the en- 
tropy of the monopole currents. According to 
Eq. JJ), fi = e 1+ f . If the monopoles make a sim- 
ple random walk on a AD hypercubic lattice then 
we get jti = 7 since there are seven choices at each 
site for the monopole current to go further. 

The entropy factor /i is shown in Figure It 
scales with b and is independent of the volume of 
the lattice. Note that in a small b region fi > 7 
because the monopole action can not be reliably 
described by the quadratic terms only |15| . 

At large b the entropy factor is smaller than 
seven. We have fitted the entropy by a function: 

^^^oo + C^-\ (3) 

where fi^, C and q are the fitting parameters. 
The fit gives ^ = 1.6(4), C = 1.7(5) and q = 
1.2(2). Fixing q = 1 in Eq. © we get ^ = 
1.15(25) and C = 2.2(1). The corresponding best 
fit curves are shown in Figure [31 

The fact that the asymptotic value of the en- 
tropy is very close to the unity in the large b limit 
may have a simple explanation. A monopole with 
a large blocking size b behaves as a classical object 
and its motion is no more a simple random walk. 
The predominant motion of the large-6 monopole 
is close to a straight line. 
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